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Spectral transfer processes in magnetohydrodynamic (MHD) turbulence are investigated 
analytically by decomposition of the velocity and magnetic fields in Fourier space into 
helical modes. Steady solutions of the dynamical system which governs the evolution of 
the helical modes are determined, and a stability analysis of these solutions is carried 
out. The interpretation of the analysis is that unstable solutions lead to energy transfer 
between the interacting modes while stable solutions do not. From this, a dependence of 
possible interscale energy and helicity transfers on the helicities of the interacting modes is 
derived. As expected from the inverse cascade of magnetic helicity in 3D MHD turbulence, 
mode interactions with like helicities lead to transfer of energy and magnetic helicity to 
smaller wavenumbers. However, some interactions of modes with unlike helicities also 
contribute to an inverse energy transfer. As such, an inverse energy cascade for nonhelical 
magnetic helds is shown to be possible. Furthermore, it is found that high values of 
the cross-helicity may have an asymmetric effect on forward and reverse transfer of 
energy, where forward transfer is more quenched in regions of high cross-helicity than 
reverse transfer. This conforms with recent observations of solar wind turbulence. For 
specific helical interactions the relation to dynamo action is established. The present 
analysis provides new theoretical insights into physical processes where inverse cascade 
and dynamo action are involved, such as the evolution of cosmological and astrophysical 
magnetic fields and laboratory plasmas. 


1. Introduction 

Since the contributions by Richardson and Kolmogorov it is well established that 
the average transfer of kinetic energy occurring in isotropic non-conducting turbulent 
fluids in three dimensions proceeds from the large scales to the small scales, or, in 
the Fourier representation, from small to large wavenumbers (see e. g. Frisch (1995); 
McComb (2014)). However, many turbulent flows occurring in nature and/or in industrial 
applications deviate from this behaviour, showing some backwards energy transfer, such 
as rotating flows (Mininni et al. 2009), two-dimensional flows e. g. (Kraichnan 1967; 
Boffetta & Musacchio 2010; Mininni & Pouquet 2013) as well as flows doted with 
polymeric additives (Dubief et al. 2013). 

Electrically conducting turbulent flows also show a variety of phenomena resulting in a 
transfer of energy from the small scales to the large scales (Biskamp 1993). One of these 
is the inverse cascade of magnetic helicity, first proposed by Frisch et al. (1975) after 
the derivation of absolute equilibrium spectra for magnetic and kinetic energies as well 
as cross and magnetic helicities and subsequently confirmed numerically (Pouquet et al. 
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1976; Pouquet & Patterson 1978). By virtue of the realisability condition, which states 
that magnetic energy at a given wavenumber k is bounded from below by k/2 times the 
magnitude of the magnetic helicity (see e. g. Moffatt (1978); Biskamp (1993)), it also 
drives a transfer of magnetic energy from the small to the large scales. Another process, 
which is of particular importance in astrophysical fluid dynamics due to the generation of 
magnetic fields of astrophysical bodies, is the large-scale dynamo, by which a magnetic 
field on scales larger than the typical scale of the largest eddies is generated from a 
magnetic seed field. The a-effect of mean-field electrodynamics is a classic example of a 
large-scale dynamo, and it relies on the presence of kinetic helicity (Krause & Radler 1980; 
Brandenburg 2001). Thus in both these cases energy is transferred in spectral space from 
large to small wavenumbers, and a lack of mirror symmetry (i.e. the presence of kinetic 
and/or magnetic helicity) facilitates these types of energy transfer. Moreover, results from 
recent numerical simulations (Brandenburg et al. 2015; Zrake 2014; Berera & Linkmann 
2014) show that some kind of inverse spectral transfer also occurs in conducting flows 
with vanishing magnetic and kinetic helicities. Recently, Stepanov et al. (2015) calculated 
energy and helicity transfer fluxes numerically from a helical shell model, where helicity 
and energy input was separated in order to test the influence of magnetic helicity on the 
turbulent dynamics. They found that the inverse cascade of magnetic helicity is connected 
to a transfer of purely magnetic energy to lower wavenumbers. By alternating the sign 
of helicity injected into the system it was shown that even if the average injection rate 
of magnetic helicity vanishes, the reverse transfer of magnetic energy induced by the 
instantaneous injection of magnetic helicity remains active. 

The effect of helicity on energy transfer and evolution in non-conducting turbulent 
fluids has received considerable attention (Moffatt 1969; Andre & Lesieur 1977; Pelz 
et al. 1986; Polifke & Shtilman 1989; Polifke 1991; Waleffe 1992; Chen et al. 2003a,b; 
Biferale et al. 2012, 2013; Biferale & Titi 2013) and has been studied in a variety of 
ways, e. g. using analytical methods, closure calculations, conventional direct numerical 
simulations (DNSs) and novel approaches in DNS. Waleffe (1992) decomposed the Fourier 
transform of the velocity field into eigenfunctions of the curl operator and derived 
evolution equations for these eigenfunctions by substitution of the decomposed field 
into the Navier-Stokes equation for incompressible flow. Since the nonzero eigenvalues 
of the corresponding eigenfunctions are related to the helicity of given velocity field 
mode, the evolution equations were further analysed in order to derive the dependence 
of the possible energy transfers on the helicities of the interacting modes. If the largest 
two wavenumbers of a given wavevector triad had helicities of opposite sign energy was 
transferred forward in wavenumber space, while a reverse transfer of energy became 
possible if the helicities were of the same sign. 

Biferale et al. (2012) investigated numerically whether this reverse spectral transfer 
caused by interactions of helical modes of the same sign occurs. By defining a projection 
operation on the nonlinear term the authors altered the Navier-Stokes equation so 
to ensure that only modes of, say, positive helicity were present in the system. The 
altered Navier-Stokes equation was subsequently solved numerically using the standard 
pseudospectral method in conjunction with small-scale forcing, and, as predicted by 
Waleffe’s analysis, it was found that kinetic energy was indeed transferred downwards 
in wavenumber space. This was the first observation of an inverse energy cascade in 
three-dimensional isotropic turbulence. 

In a subsequent paper (Biferale et al. 2013), the same authors forced the system at the 
large scales in order to study the predicted forward cascade of kinetic helicity, which was 
indeed observed in the simulations. Since the subset of positively helical modes does not 
transfer energy to the small scales, it was expected that the resulting dynamical system 
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would not show finite dissipation in the limit of infinite Reynolds number. Hence the 
projected Navier-Stokes equation which governs the evolution should be globally regular, 
which was subsequently proven by Biferale & Titi (2013). 

Thus, in summary, the decomposition of the Fourier transform of the turbulent velocity 
field fluctuations into helical modes has been proven to be very useful in terms of under¬ 
standing some fundamental features of turbulent flows which go beyond the established 
Kolmogorov-Richardson (direct) cascade of kinetic energy. In view of the effects of kinetic 
and magnetic helicities on the direction of energy transfer in MHD turbulence and 
inspired by the successes of the helical decomposition used in hydrodynamics, in this 
paper we propose to use the decomposition of both the magnetic and velocity fields into 
helical modes in order to perhaps shed some more light on why magnetohydrodynamic 
(MHD) turbulence shows much more transfer from the small scales to the large scales 
than turbulence in non-conducting fluids. This analysis provides deeper insight into 
fundamental transfer processes in MHD turbulence. It will help in the further theoretical 
understanding of various physical processes involving inverse cascade and dynamo action, 
such as the evolution and generation of cosmological and galactic magnetic fields. 

Before outlining the structure of this paper we pause briefly to discuss the terminology 
used. As the precise meaning of the term varies in the literature, it is not always evident 
what is meant by an inverse cascade. In the astrophysical literature, transfer of energy and 
helicity from higher to lower wavenumbers is often described as an inverse cascade (Son 
1999; Christensson et al. 2001; Cho 2011), while the fluid dynamics literature requires 
any cascade process to possess a wavenumber-independent flux (McComb 2014; Biferale 
et al. 2012; Biferale & Titi 2013; Muller et al. 2012; Alexakis et al. 2006; Brandenburg 
et al. 2015). It is thus of interest to not only classify the different types of reverse 
transfer that can occur in MHD turbulence, but also to perhaps clarify the terminology. 
Therefore we propose in this paper the general term reverse (or inverse) spectral transfer, 
which encompasses all the phenomena described above as subcategories according to their 
properties. We define it as any process that produces an increase in a spectral quantity 
(total energy, magnetic helicity, etc.) at low wavenumbers due to transfer of that quantity 
away from higher wavenumbers into smaller wavenumbers. In this framework an inverse 
cascade is a reverse spectral transfer showing constant flux of the cascading quantity 
over a certain wavenumber range. We point out that in MHD concerns have been raised 
over the use of the term ‘cascade’ (Muller et al. 2012), as it may be understood to imply 
energy (or magnetic helicity) transfer mainly due to local interactions, which might not 
be the case in MHD turbulence (Alexakis et al. 2005; Debliquy et al. 2005; Cho 2010; 
Muller et al. 2012). We will come back to this point in the discussion section of this 
paper. 

This paper is organised as follows. First we explain the helical decomposition of the 
velocity and magnetic fields in sec. 2 and use this decomposition to outline the derivation 
of the evolution equations for the helical coefficients from the MHD equations in sec. 3 
following the work by Lessinnes et al. (2009). This leads to a system of coupled ordinary 
differential equations (ODEs) describing the interaction of modes corresponding to a 
single triad of wavenumbers. In sections 4 and 5 the linear stability of steady solutions of 
this system is examined. Section 6 explains the assumptions used to interpret results from 
the stability analysis followed by a presentation of Hrst results on energy transfers, which 
are used in sec. 7 to calculate the contribution to inertial range energy and magnetic 
helicity cascades. In sec. 8 we analyse specihc cases where the MHD equations can be 
decoupled and relate our analysis to large- and small-scale kinematic dynamo results 
such as the a-effect (Moffatt 1978; Parker 1979; Krause & Radler 1980). Our results 
are summarised and discussed in sec. 9 in the context of numerical simulations of MHD 
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turbulence and observations of turbulence in the solar wind, at this point we also provide 
suggestions for further work. 


2. Problem statement and basic equations 

In this paper we will be studying the interscale energy and helicity transfers that 
govern the dynamics of MHD turbulence in the absence of a background magnetic field. 
We consider the turbulent flow to be incompressible and make no explicit assumptions 
about the value of the magnetic Prandtl number. The MHD equations for incompressible 
flow are 

dtu = - VP — (u ■ V)u H —(V X 5) X b + lyAu , (2-1) 

P P 

dtb = (6 • V)m — (u ■ V)b + r/Ab , (2-2) 

V • a = 0 and V • 6 = 0 , (2.3) 


where u denotes the velocity field, b the magnetic induction expressed in Alfven units, 
ly the kinematic viscosity, y the resistivity, P the pressure and p the density, which is set 
to unity for convenience. 

For simplicity at first we consider periodic boundary conditions on a domain 17 = 
[0,P]^ C thus working with the discrete Fourier transformed MHD equations 


V P 


{dt + iyk‘^)u{k) = - FT 


+ ^ -{ip X u{p)) X u{q) + {ip X b{p)) X b{q) 

k+p+q=0 

{dt + vk'^)b{k) = ik X ^ u{p) x b{q) , 

0=k+p+q 


(2.4) 

(2.5) 


where FT denotes the three-dimensional Fourier transform as a linear operator acting 
on L^{[2) functions, u and b denote the Fourier transforms of the velocity and magnetic 
fields respectively and the inertial term {u ■ V)it in the momentum equation (2.1) has 
been written in rotational form {u ■ V)m = (V x it) x tt -|- 0.5V|ttp. In order to determine 
the contribution of specific interactions to the fluxes of magnetic helicity and magnetic 
energy, eventually we will formally take the limit P —)■ oo in sec. 7, necessarily assuming 
that the relevant functions are then well-behaved at infinity to ensure the convergence 
of the respective Fourier integrals f. 


2.1. Helical decomposition 

The decomposition of the Fourier transform of a solenoidal vector field in circularly 
polarised waves as proposed by Constantin & Majda (1988) has been used in several 
investigations of hydrodynamic turbulence (Waleffe 1992; Biferale et al. 2012; Biferale 
& Titi 2013) in order to establish the properties of energy transfer depending on the 
kinetic helicity. For conciseness we only review the fundamental properties of the helical 
decomposition and refer to the relevant literature for details and derivations. 

The action of the curl operator on square integrable real vector field v{x) can be viewed 
in spectral space as the action of a linear operator acting on the Fourier transform v{k) 

f A discussion of this point can be found in the book by Titchmarsh (1939) 
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4 : ^ 

v{k) —ik X v{k) . 

As such the linear operator /fe(-) = ik x (•) has a set of linearly independent eigenvectors 
defining a basis of C^, thus v{k) can be expanded in this basis. That is, it can be expressed 
as a linear combination of eigenvectors ik, h^ik) and h^ik) of the curl operator 4, where 


%k X — Sf^hhg^ , 

(2.6) 

-ik X = Skkhl^ , 

(2.7) 


Sfc = ±1 and Skk = ±fc are the nonzero eigenvalues of the curl operator in spectral space 
t, and * denotes the complex conjugate. The complex eigenvectors are orthogonal to each 
other and are fully helical. They are normalised to unit vectors for the remainder of this 
paper. 

Since u{k) and b{k) are solenoidal, they can be expressed in terms of only 

u{k,t) = u-{k,t)h^{k) + u+{k,t)h+{k) = '^Us^{k,t)hs^{k) , (2.8) 

Sk 


h{k,t) = b-{k,t)h^{k)+ b+{k,t)h+{k) ='^bs,^{k,t)hsk{k) , (2.9) 

where Ua^, and bg^ are complex coefficients. These coefficients can be calculated by taking 
the inner product of the basis vectors with the appropriate fields 


and 




bs^{k,t) = 


■ u{k,t) 

hg^k) ■ hl^{k) ’ 

(2.10) 

K^ik)-bik,t) 
hg,{k) ■ hl^ik) ■ 

(2.11) 


In order to keep the derivations concise we suppress the dependence of the coefficients 
on time and wavevector from now on in the notation. 

The helical decompositon of a solenoidal vector field was first applied to incompressible 
MHD flows by Lessinnes et al. (2009), who derived a dynamical system in Fourier space 
describing helical triadic interactions in MHD. This system was subsequently used to 
construct a helical shell model of MHD turbulence. In the following section we briefly 
review the derivation carried out by Lessinnes et al. (2009). 


3. The evolution of the helical modes 

The equations describing the evolution of the helical coefficients and bg^, are 
derived by substituting the decompositions (2.8) and (2.9) into the MHD equations for 
incompressible flow and then taking the inner product with hg^, on both sides of the 

f Note that the curl operator can have eigenvectors with nonzero eigenvalues, as it involves the 
cross product of two complex vectors. This is not necessarily orthogonal to the plane spanned by 
the two complex vectors, instead it is orthogonal to the plane spanned by the complex conjugates 
of the two vectors. 
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respective equations. The resulting evolution equation for the helical coefficient Ug,. is 


{dt + vk'^)ug, = ^ 


-FT 


¥ E 


V ( 

— {ip X u{p)) X u{q) + {ip x b{p)) x b{q) 


k+p+q=0 


= “ 9 E E ~ 


Sp,Sq 0—k-^p-^q 


hi xh*- hi 


(3.1) 


where the dummy variables p and q were exchanged in order to symmetrise the mo¬ 
mentum equation with respect to p and q and thus to obtain the factor {spp — Sqq)/2. 
Following an analogous procedure (Lessinnes et al. 2009) for the helical coefficient bgf. of 
the magnetic field leads to 


h 


{dt + qe)bg, = '-^ 


ik X E u{p) X b{q) 

0=fc+p+q 


_ Skk 


E 

Sp,Sq 0 — k-\-p-\-Q 


hi X hi 


hi 


{ul bl -bl ul) 


(3.2) 


In order to study the interaction of helical modes, that is the evolution of the helical 
coefficients due to the mode coupling only, the diffusivities are from now on omitted. 
For a given triad k,p, q of wavevectors, expressions for the first time-derivatives of each 
helical coefficient are obtained from (3.1) and (3.2) and from the corresponding equations 
for bgp , , Usj, and Ug^ ■ This leads to the following system of coupled ODEs describing 

the evolution of the helical coefficients in a single triad interaction 


dtUg^ = {spP - Sgq) gkpq {u*g^u*g^ - bl^blj , 

dtUg^ = {sqq - Skk) gkpq {K^ul^ - blblj , 

dtUg^ = {skk - Spp) gkpq {Ug^ul^ - bl^blj , (3.3) 

dtbg^ = -Skk gkpq {u*g^bl^ - bl^ul^) , 

dtbg^ = -Sppgkpq {u*g^bl^ - biy^J , 

dtbg^ = -Sqq gkpq {ul^bl^ - bl^ulj , (3.4) 

where the geometric factor 

9kpq =^ Kf K, , ( 3 . 5 ) 

is introduced for conciseness, following Waleffe (1992) and Lessinnes et al. (2009). It can 
also be written as 

gkpq = ^(Sfefc + SpP + Sqq) , (3.6) 

where a is a wavenumber-dependent real number determined by the orientation of the 
triad and N a factor depending on the shape of the triad. Further details and a derivation 
of (3.6) can be found in the paper by Waleffe (1992). 

The three ideal invariants total energy, magnetic helicity and cross-helicity are dehned 
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respectively as 

Etot ^ E M ’ (3.7) 

k k,Sk 

H^ag=J2(^{k)bi-k) = J2"^\bsj\ (3.8) 

k k.Sk 

H, =J2{Mk)b{-k)) = J2^e {u.Xj , (3.9) 

k k,Sk 

where a denotes the vector potential of the magnetic field, Re the real part of a complex 
number and angle brackets an ensemble average. They are conserved in single triad 
interactions (Lessinnes et al. 2009). 


4. Stability of steady solutions 

Examining the linear stability of steady solutions of the system (3.3)-(3.4) can reveal 
the influence which the helicities of the interacting modes have on the interscale transfer 
of a given quantity of interest. 

The system (3.3) without a magnetic field (that is for bs = 0) was analysed by 
Waleffe (1992) with respect to the linear stability of its steady solutions. Linearly 
unstable solutions were found depending on the helicities of the interacting modes. 
This result was then interpreted following the instability assumption inspired by the 
formal analogy to rigid-body rotation, where rotation around the axis of middle inertia 
is unstable. The existence of a linearly unstable solution involving a velocity field mode 
u is interpreted as the M-mode losing energy to the other two modes it interacts with. An 
equivalent assumption had already been used by Kraichnan (1967) for two-dimensional 
hydrodynamic turbulence. In the remainder of this paper we take a similar approach 
and investigate the linear stability of steady solutions of the system (3.3)-(3.4) in view 
of possible applications to spectral transfer processes in MHD and in particular for the 
inverse transfers of total energy and magnetic helicity. In principle, a similar analysis 
could be carried out for the remaining ideal invariant, the cross-helicity. 

4.1. The steady solutions 

The system (3.3)-(3.4) of six coupled DDEs has several steady solutions one can 
linearise about. To simplify the notation, a (formal) solution of the system (3.3)-(3.4) 
consisting of helical u- and 6-field modes interacting in a given triad k, p, q is written as: 

iusk,Us^,Us,]bsk,bs^,bsJ . (4.1) 

In order to find the steady solutions of the system (3.3)-(3.4), we assume (without 
loss of generality) that the middle components and Ug^ = Ug^ are constant in 

time. Then (3.3) and (3.4) require the other four components to vanish by the following 
argument. A steady solution requires dfUg^, = 0, and the only way that this can happen 
nontrivially is if both products u*^u*^ and 6*^6*^ vanish f. Since Ug^ = Ug^ is constant 
in time, Ug^ = 0 and similarly bg^ = 0. This leaves us with 

{usk ^ E^sp j 0; : Bg ^, 0) . 

f This requires assuming that no cancellations occur. However, the occurrence of cancellations 
would require the system to be in a specific state, which is unlikely to happen frequently in a 
chaotic system. 
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Applying the same argument to dtUs^, it follows that Us^. and bs^. must also vanish. 
Therefore a steady solution of the system (3.3)-(3.4) has the form 

It can now be checked for consistency that dtbg = 0 for k, p and q as well. Therefore 
the solution is steady for the magnetic held and for the velocity held alike. Asides 
from the just explained example, steady solutions of the form {Us^,0,0; Bs^,0,0) and 
(0,0, Usg ; 0,0, Bs^) are obtained in the same way. 

Thus the steady solutions of (3.3)-(3.4) are of the same form as for the hydrodynamic 
case (Waleffe 1992), where at least two of the three interacting modes vanish. However, 
there are two special cases: one where the magnetic held component Bg also vanishes, 
while Us 0 and the other, where the velocity held component Ug vanishes, while Bg ^ 0. 
The former case may perhaps be connected to a dynamo process. At this point we note 
that for the kinematic dynamo, where the back-reaction of the magnetic held on the 
velocity held can be neglected, the (linear) stability of the velocity held coefficients Ug is 
only determined by hydrodynamic interactions. We will come back to this point in sec. 8. 


4.2. Linear stability analysis 

In order to assess whether a given steady solution is linearly stable in our particular 
setting, we assume without loss of generality that the coefficients Ug^ and bg^ correspond¬ 
ing to wavevector p are nonzero and constant in time, that is, we study the linear stability 
of the solution (0, Us^,0] 0, Bg^,0) with respect to inhnitesimal perturbations of the four 
modes that had been set to zero. As the hrst-order equations involve the coupling of 
all three modes of a given triad, little information can be obtained from them at hrst 
sight. Therefore we pass to second-order time-derivatives of the evolution equations for 
the perturbations Ug^, bg^, Ug^ and bg^. Taking time-derivatives on both sides of (3.3)- 
(3.4) and subsequently substituting any occurrence of a hrst-order time-derivative on the 
right-hand side by the appropriate evolution equation, we obtain 


d^bg. 


\gkpq P(SpP- Sgg) 
\9kpq I {SpP Sqq) 


[((Sfefc - Spp)\UgJ^ + Sqq \BgJ^)ug^] 
[((sfcA: - Spp)U:^Bg^ + Sqq Ug^B*Jbg,] , 


\9kpq \^Skk {sqq Ug^Bg^ -f {skk 

\9kpq I Skk I^SpI 4” {^kk 


- Spp)Ug^Bl^)ug^ 
SpP)\Bs^?)bg^] , 


(4.2) 


(4.3) 


where terms of second order in small quantities (such as e.g. uj.) have been neglected. 
Note that these equation do not depend on modes at wavenumber q. The evolution 
equations of the helical coefficients Mg, and bg^ can be obtained similarly and show no 
dependence on k, therefore we restrict our attention to the evolution of Ug^. and hg^. 
The system (4.2) and (4.3) can be written as a matrix ODE 


where x = {ug ^, bg ^) and the matrix elements are 

a = \9kpq \^{spP- Sqq)[{skk - Spp)\Ug^\^ + Sqq \Bg^\^] , 

P = -\9kpq ?{spP - Sqq) [{skk - Spp)Ug^Bg^ + Sqq Ug^B*gJ , 
7 = \9kpq fskk Ug^Bg^ + {skk - Spp)Ug^BgJ , 

^ ~ \9kpq I ^kk \_^qQ I^SpI 4“ {^kk ^pP)\Bsp\ ] . 


(4.4) 


(4.5) 

(4.6) 

(4.7) 

(4.8) 
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The linear stability of this system can be determined from the eigenvalues Ai and A 2 of 
the matrix in (4.4). These eigenvalues depend not only on the helicities of the interacting 
modes and on the magnitudes of Us^ and Bs^ relative to each other, but also on the 
alignment between the magnetic and velocity field modes at wavevector p, that is, on 
the cross-helicity. For a given steady solution to be unstable the perturbations have to 
be exponentially growing, and so at least one of the eigenvalues (for i = 1,2) must 
have a positive real part. We will now assess under which conditions this is possible. 
The eigenvalues Ai {i = 1,2) are given by 

, a + 5 I{ct + S)"^ ^ 

Ai ,2 = ^—± y- ^ - a6 + p-f. (4.9) 

For convenience define 

X = —and Q = a6 — Pj , (4-10) 

such that 

X = - \UsJ'^[skksqq + {skk - Spp){sqq - Spp)] 

- \Bsj'^[Skk{Skk- Spp) + SqqiSqQ- Spp)] , (4.11) 

and 

Q = \9kpq\^Skksqq{skk - Spp){sqq - Spp) + 2\UsJ‘^\BsJ'^ - 4i/c(p)^)) , 

(4.12) 

hence the cross-helicity H^p) enters the dynamics through the parameter Q. The 
derivation of (4.12) can be found in appendix A. Note that the term + \Bsp\^ + 

2|17sj,P|i3spP — 4i7c(p)^ is always positive, regardless of the value of He since \Hc{p)\ < 
\Usj,\\BsJ, thus the sign of Q is determined by the helicities of the interacting modes and 
the wavenumber ordering. 

The eigenvalues Ai can now be written more concisely as 

An 2 = x± \/x'^ - Q , (4-13) 

therefore the possibility of finding exponential solutions of the system (4.4) depends on 
the values of x and Q. Apart from the trivial case, where x = 0 and Q = 0, there is 
only one case for which no linear instability occurs: this is if x < 0 and |x| > \\/x^ — Q\, 
since then -y/Ai and are imaginary numbers allowing only oscillatory solutions of the 
matrix ODE (4.4). All other cases lead to exponentially growing as well as exponentially 
decaying solutions. 

Cases in which x > 0 and Q < 0 result in the largest eigenvalues and thus in the fastest 
growing exponential solution. These cases are therefore of special interest, as within the 
framework of the instability assumption they may lead to the largest energy transfer and 
thus can yield information about which combination of parameters facilitates most of the 
energy and helicity transfers. We will consider this point in further detail in sec. 7. 

As can be seen from the structure of the terms x and Q, the relative magnitudes 
and the ordering of the wavenumbers in a given triad will influence the stability of 
steady solutions. In view of the continuous interest in nonlocality of interactions in MHD 
turbulence (Brandenburg 2001; Alexakis et al. 2005; Debliquy et al. 2005; Cho 2010; 
Muller et al. 2012), we point towards specific results for local and nonlocal interactions 
where appropriate. Following Waleffe (1992), for wavenumbers ordered k < p < q, the 
nonlocal limit is defined as fc << p ~ < 7 , while local interactions are characterised by 

k p q. 
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5. Instability and helical interactions 

Since s = ±1, interactions between helical modes which all have helicities of opposite 
signs are not possible, and at least two modes will always have helicities of the same sign. 
Therefore we have four classes of possible helicity combinations 

Sfc = Sp ^ Sq , Sk — ^ ^ ~ ^p and Sfc — Sq = Sp , 

each of which occur twice as s can take the values ±1. These four possible (classes of) 
combinations are now studied on a case-by-case approach in order to determine when a 
certain combination of helicities leads to exponentially growing solutions of the system 
(4.4). 

5.1. The case Sk = Sq ^ Sp 

Since the expressions in square brackets of (4.11) become 

kq + {k + p){q + p) > 0 and k{k + p) + q{q + p) > 0 , (5.1) 

one obtains x = {a + 5)/2 < 0. For an unstable solution \x\ < \yj x"^ — Q\, however, we 
obtain Q > 0 since 

Q ~ SkSqkq{skk - Spp){sqq - Spp) , (5.2) 

which is positive for Sk = Sq ^ Sp. Furthermore we obtain Q < x^ (see appendix B) 
and thus \x\ > l^/x"^ — Q\, which results in imaginary eigenvalues of the matrix in (4.4). 
Therefore we do not obtain unstable solutions for the case Sk = Sq ^ Sp, and this is 
independent of the ordering of the wavenumbers k,p and q. Note that this implies that 
exponentially growing solutions of (4.4) are impossible if the perturbations , Us ^, 
and hs^ have helicities opposite to the helicities of the modes Us^ and constituting 
the steady solution. 

For the remaining helicity combinations, which do result in unstable solutions, the 
ordering of wavenumbers matters. The arguments used to decide whether or not an 
exponentially growing solution becomes possible are similar to the procedure employed 
for the case Sk = Sq ^ Sp described above. 

5.2. The case Sk ^ Sp = Sq 

In this case we obtain 

Qkq{k + p){q - p) (5.3) 

and 

^ = -\El^^\UgJ^[-kq - {k + p){q-p)] - \BgJ^[k{k + p) + q{q - p)] . (5.4) 

The stability of a steady solution depends on the signs of these terms which in turn 
depend on wavenumber ordering, cross-helicity and the ratio |C^sp |- 

• For k < p < q we obtain unstable solutions if \UsJ > |i3sp|, since then x > 0. 
For |i?sp| > |C/sp| unstable solutions are still possible, provided Hc{p) is small and \BsJ 
not much larger than |I7sp|- Thus in regions of large cross-helicity unstable solutions 
only occur for weak magnetic fields. The method by which these resnlts are obtained is 
explained in appendix C. 

For nonlocal interactions (fc << p ~ g) we obtain Q = 0 and the sign of x determines 
whether unstable solutions occur. The term x is now of the form 

''-^^kq{\Ugf -\Bgf) , 


X ~ 


(5.5) 
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hence nonlocal interactions lead to unstable solutions if \BsJ < \UsJ. 

• For k < q < p, Q will become negative, leading to unstable solutions regardless of 
the ratio \UsJ/\BsJ and the value of Hc{p). 

• For p < k < q again we obtain unstable solutions if \Usj,\ > \BsJ, since then 
a: > 0. For \BsJ > \Usj,\ unstable solutions are still possible, provided Hc{p) is small 
and \Bsp\/\Usp\ not >> 1 (see appendix C). Nonlocal interactions {p « k ~ q) lead to 
unstable solutions if \Us^ \ > iBgJ, because then 

x^\gkp,\^k^{\Usf-\B^f)>0. (5.6) 

In summary, a given steady solution in this case is more likely to be stable if the nonzero 
mode is at medium or low wavenumbers in regions of high cross-helicity. 

5.3. The case Sk = Sp = Sq 

In this case we obtain 

Qkq{k - p){q - p) (5.7) 

and 

x = -'^-^^^^^\Us^\’^[kq + {k-p){q-p)\- \Bs^’^[k{k - p) + q{q - p)] . (5.8) 

• For k < p < q we obtain Q < 0 and thus x + \/x^ — Q ^ 0, leading to exponentially 
growing solutions independent of i7c(p) and the ratio \Usp\/\Bs^\. We note that both 
velocity and magnetic field modes have positive and negative contributions to the sign 
of X. This is of interest since if a: > 0 the resulting eigenvalue would be larger and thus 
the solution would grow faster. However, in this case this cannot be determined from the 
ratio \Usp\/\Bg^\ and thus there is little information about what contributes to a faster 
growing exponential and thus to a more unstable solution. 

For both local {k ~ p ~ q) and nonlocal {k « p ~ q) interactions we obtain Q = 0 and 
the sign of x determines whether unstable solutions occur. For the nonlocal case only the 
magnetic field term is positive, and x has the form 

x^'^-^^kq{\BgJ^-\Ugf) . (5.9) 

leading to unstable solutions if |Ssp| > iC^spI) while for local interactions no instability 
occurs as the only term in x that does not vanish is —\gkpq\^\Usp\^kq < 0. 

• For k < q < p, the possibility of exponentially growing solutions depends on the ratio 
\Usj,\/\Bsj, \ and on the relative magnitudes of the wavenumbers k,p and q, as now Q > 0. 
Since the magnetic field term in x is now positive, instabilities occur for iH^p |/|i?sp | < 1. 
If IIfsp I/1-Ssp I > 1 it depends also on the cross-helicity whether instabilities occur. For 
maximal Hf.{p) we obtain — Q > 0, hence the solutions will be stable. If Hc{p) = 0 and 
\Us \/\Bs I is not too small, instabilities will occur, depending also on the shape of the 
triad (see appendix C for further details). In general, the smaller \Us \/\Bs \ the more 
unstable is the solution. 

• For p < k < q we obtain a; < 0 and Q > 0, furthermore a;^ — Q > 0 independent 
of |Ifsp|/|I?sp| and Hc{p) (see appendix C), thus no unstable solutions occur. Nonlocal 
interactions {p « k ~ q) do not lead to unstable solutions, since 


x^-\gupqnk^-kp]{\UgX + \BgX)<Q. 


( 5 . 10 ) 
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Helicities | He \ constraint | stability 


Sk ^ Sq ~ Sp 

n/a 

\Usp\ > \Bsp\ 

unstable 


max 

IB.pl > \Usp\ 

stable 


0 

\Bsp\/\Usp\ not » 1 

unstable 

Sk = Sp ^ Sq \ 

n/a 

n/a 

unstable 


Sfc = Sg = Sp I n/a I n/a | unstable 


Table 1: Summary of possible unstable solutions for the middle wavenumber modes k < 
p<q. 


5.4. The case Sk = Sp ^ Sq 

The terms determining the stability in this case are 

Qkq{k - p){q + p) (5.11) 

and 

X = \u^j'^[-kq- {k-p)iq+p)] - \B,j'^[k{k-p) + q{q+p)] . (5.12) 

• For k < p < q unstable solutions occur independent of the ratio \UsJ/\Bsj,\, and 
since both magnetic and velocity field terms have positive and negative contributions to 
the sign of x, we are in a similar situation to the previous case. However, in the present 
case Q ~ 0 only for local (fc ~ p ~ g) interactions. It is now the velocity field term 
\ 9 kpq\‘^\Usp\'^kq > 0 which ensures that exponentially growing solutions exist for local 
interactions provided \Usp \ > 2\Bsp\. 

• For k < q < p the result is the same, since reversing the relative ordering of p and q 
does not change the sign of Q. That is, exponentially growing solutions occur. 

• For p < k < q the term Q is positive and the term proportional to is positive 

while the term proportional to is negative. Thus instabilities occur if > 

1. For \Usp\/\Bsp\ < 1, the occurrence of instabilities depends on the value of Hc{p). If 
Hc{p) is maximal and magnetic and velocity field are fully aligned, then the solutions 
are stable. For zero cross-helicity and \Bs \ being not much larger than \Us |, solutions 
are unstable (see appendix C). 

This type of helicity combination is another possibility for nonlocal interactions of the 
type p << A: ~ g leading to exponentially growing solutions if \Usp \ > |A?sp|, since then 

x^\gkpq\^k^{\Usp\^-\B^pn>0. (5.13) 

The results of the dependence of the occurrence of unstable solutions on combinations 
of helicities, wavenumber ordering, relative magnitudes of the u and b modes and cross- 
helicities at wavenumber p are summarised in tables 1-3. 


6. Energy transfers and the instability assumption 

In order to use the results of the previous section to derive results for the transfers 
of the ideal invariants total energy Etot and magnetic helicity Hmag, we invoke the 
instability assumption (Waleffe 1992). Generalised to MHD turbulence, this assumption 
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Helicities 

Ho 

constraint 

stability 

Sk ^ Sq = Sp 

n/a 

n/a 

unstable 

$k — Sp Sq 

n/a 

n/a 

unstable 

Sk — Sq = Sp 

n/a 

IB.pl > \Us,\ 

unstable 


max 

\UsJ > \BgJ 

stable 


0 

\UsJ/\BsJ not » 1 

unstable 

Table 2: Summary of possible unstable solutions for the largest wavenumber modes k < 
q<p. 


Helicities 

Ho 

constraint 

stability 

Sk ^ Sq — Sp 

n/a 

\Us,\ > \BsJ 

unstable 


max 

|B.J > \UsJ 

stable 


0 

\BsJ/\Usp\ not » 1 

unstable 

Sk ~ Sp Sq 

n/a 

\UsJ > \BsJ 

unstable 


max 

|B.J > \UsJ 

stable 


0 

\BsJ/\UsJ not » 1 

unstable 

Sk ~ Sq = Sp 

n/a 

n/a 

stable 

Table 3: Summary of possible unstable solutions for the smallest wavenumber modes 
p < k < q. 


asserts that energy is transferred away from modes whose evolution equation for the 
helicity coefficient is linearly unstable, into the other two modes it is coupled to by a 
triad interaction given through the system (3.3)-(3.4). 

Therefore the results of the stability analysis determine whether a given helicity 
combination mainly contributes to forward or reverse transfer of energy. That is, if a 
steady solution at wavenumber p is unstable and energy is transferred away from 
and Usj, into the modes they interact with (note that Bg^ and Ug^ do not interact with 
each other directly), then the wavenumber ordering k < q < p results in reverse transfer 
of energy, while p < k < q results in forward transfer and k < p < q in a split transfer 
with contributions to forward and reverse directions of energy transfer. 

Several immediate results can be deduced from the summary of the stability analysis 
for the different helicity combinations presented in tables 1-3. First, unlike in non¬ 
conducting fluids modes corresponding to the largest wavenumber in a given triad can 
be unstable, leading to more possibilites for reverse spectral energy transfer in MHD 
compared to hydrodynamics. Second, all three helicities influence the direction of energy 
transfers, and reverse transfers are also possible for cases of unlike helicities. Third, 
forward transfers appear to be more quenched in regions of high cross-helicity than 
reverse transfers. Fourth, very nonlocal triads contribute mainly to reverse transfers in 
magnetically dominated systems through interactions of modes with like helicity. They 
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only contribute to forward transfers through interactions of modes with unlike helicity 
and mostly if the kinetic energy is larger than the magnetic energy. 

Therefore we obtain that reverse spectral transfer becomes much more likely in MHD 
turbulence than in turbulence of non-conducting fluids, which reflects the predictions 
from absolute equilibrium spectra (Frisch et al. 1975; Zhu et al. 2013) and the well- 
established numerical results on inverse cascades, and more generally reverse transfer, in 
MHD turbulence (Pouquet et al. 1976; Pouquet & Patterson 1978; Balsara & Pouquet 
1999; Alexakis et al. 2006; Brandenburg 2001; Muller et al. 2012; Berera & Linkmann 
2014; Brandenburg et al. 2015). 

We note that the transfer directions deduced so far may or may not contribute to 
forward and inverse cascades of energy and magnetic helicity, as no information on the 
constancy, or otherwise, of the fluxes of these quantities through a given wavenumber is 
available at this point. The aim of the next section is to determine the contribution of 
the individual transfers to energy and magnetic helicity cascades. 


7. Transfer and cascades of total energy and magnetic helicity 

In order to determine the contribution of a given interaction of helical modes to 
energy and magnetic helicity cascades, the fluxes of these quantities need to be calculated 
and studied in the respective inertial ranges where they are wavenumber-independent. 
However, several technical details need to be discussed before we can proceed to this 
calculation. 

In the discrete Fourier representation the evolution equations of the kinetic and 
magnetic energy spectra Ekin{k) and Emag{k) are obtained by multiplying the relevant 
equations in the system (3.3) by and b*^, respectively, then summing over all triads 
and helicity combinations and finally carrying out shell- and ensemble averages. For the 
kinetic energy spectrum this leads to 

.. zl 8 

dtEkin{k') = — EE - (7-1) 

p,q i—1 

where denotes a sum over all wavenumbers p and q whose wavevectors p and q 

form a triad with k such that k + p + q = 0 and the superscript (i) labels the eight 
possible helicity combinations. The transfer terms in this equation are given by 

= (Spp- Sqq) gkpqius^Us^Us^) +C.C. , (7.2) 

Sik,p,q) 

and 

t^'Lpik^P.q) =-{spp-Sqq) Y 9 kpq{uskBsj,hs^) + c.c. , (7.3) 

S(k,p,q) 

where S{k,p, q) indicates a summation over all wavevectors in shells of radius k,p and q 
and c.c. denotes the complex conjugate. For the magnetic energy spectrum one obtains 

, Zi 8 

dtE^ag (fc) = 2 E E ik,P,q) , (7-4) 

p,q i—1 

where 

^magik:P: q) = ~Skk ^ ) QkpqibskBg^Us^ — bg^Us^bs^) -|- C.C. . (7-5) 

S{k,p,q) 
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The evolution equation for the total energy spectrum E{k) = Ekin{k) + Emag{k) is 
given by the sum of the respective evolution equations for Ef^in(k) and E^agik) 

, (7.6) 

and total energy transfer term t^'^'>{k,p,q) therefore consists of three types of transfers 

t^"Hk,p,q) = t%^jy{k,p,q) + t^l^p{k,p,q) + tl^^g{k,p,q) . (7.7) 

The term t^^jj{k,p, q) denotes purely hydrodynamic transfer due to the coupling of the 
velocity field to itself, t^l\p{k,p,q) the contribution due to the Lorentz force acting on 
the fluid and tmag{k,p, q) the contributions due to advection of the magnetic field by the 

flow and conversion of kinetic to magnetic energy, that is, due to dynamo action. In real 

space the nonlinear term V x (m x 6 ) corresponding to the magnetic transfer term can be 
split into an advective term (tt • V)b and a dynamo term [b ■ V)t 6 , however this splitting 
is obscured in Fourier space. 

These terms are still written in the discrete Fourier representation of the magnetic 
and velocity fields. However, the calculation of the energy and magnetic helicity fluxes 
requires a continuous Fourier representation. The continuous transfer terms are given in 
terms of Fourier integrals and can formally be obtained by taking the period L to infinity, 
assuming that the respective integrals are well-defined. The sums then become integrals 
and the continuous counterpart of e.g. the hydrodynamic transfer term t^^j^ becomes 

THD(k,P,q)dk dp dq = lim t^^jg{k,p,q) 

= {spp-Sqq)l dk dp dq gkpq{us^Us^Us^)+C.C. . (7.8) 

J\k\=k d\p\=p d\q\=q 

The transfer terms and Tmag are defined analogously. 


7.1. Total energy transfer 

In the absence of dissipation the total energy is conserved and the transfer term 
T{k,p,q) in the spectral evolution equation of the total energy redistributes energy 
between the Fourier modes and vanishes if integrated over all space. Therefore the flux 
of total energy through wavenumber k due to a given interaction (i), 

pk poo poo 

i7W(fc) = -/ dk' / T^'\k',p,q)dpdq , (7.9) 

^0 J k J k 

can be written as the sum of two contributions: the flux of total energy into all modes 
at wavenumber k' due to triads with p, q < k < k' minus the flux of total energy into all 
modes at k' due to triads with k' < k < p,q 


-| poo pk pk -1 pk poo poc) 

n^^{k) = - dk' / T^'\k',p,q)dpdq-- dk' / T^'\k',p,q)dpdq . 
^ Jk Jo Jo ^ Jo Jk Jk 


(7.10) 

We now follow the procedure introduced by Waleffe (1992) in order to render the two 
integrals in (7.10) independent of k. This is achieved using a scaling argument, where the 
two integrals are treated separately. For conciseness we only outline the procedure briefly 
for the first integral on the RHS of (7.10) and refer to the original work of Waleffe (1992) 
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for the full derivation. The aim is to express the transfer function in the first integral on 
the RHS of (7.10) in terms of new variables 


q k' k 

V = -, w = —, u = — , 

p p p 


(7.11) 


in order to remove k from the integration limits. Since ,p, q) may scale differently 

compared to T^p{k',p,q) and Tmag{k',p,q), the term T^^'> {k',p, q) in (7.10) must be re¬ 
placed by the individual transfer terms. The transfer terms are now expressed individually 
in terms of the new variables u, v and w 

T%{k',p,q) = p~I^T^h^jj{w,1,v) = r4*^(w,l,u) , (7.12) 


Ti%ik',p,q) = p = 


Tlf{w, 1,w) , 


(7.13) 


and 


Tj;,l{k',p,q)=p-^'T(^l{w,l,v) = , (7.14) 


where f3 is related to the exponent of the kinetic energy spectrum provided it has a 
power-law dependence on fc, while the exponent jd' is related to the exponents of the 
kinetic and magnetic energy spectra as explained in further detail in appendix D. The 
first term on the RHS of (7.10) then becomes 


1 

2 


n k 

T^'^\k',p, q) dp dq 


= k^-P / dv 


r*l+i? 


dw I du \ — 
. u 


4-/3 




U3-/: 




dv 


dw I du \ — 
. u 


4-^' . 


1 , v) + T^lJw, 1 , v) 


(7.15) 


The second term on the RHS of (7.10) can be treated similarly (Waleffe 1992), and we 
obtain 



(7.16) 


Combining the two results and integrating over u leads to the following expression for 
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the total energy transfer flux 



(7.17) 


where O^w^l^w^l + f due to the triad geometry. This now enables us to 
study the contribution to the total energy transfer from a given interaction (i), where 
the scaling of the magnetic and kinetic energy spectra will influence the transfer through 
the exponents P and P'. In the inertial range of total energy the energy transfer flux 
through a given wavenumber k does not depend on that wavenumber, which leads to 
the characteristic values of the scaling exponents P' = P = 3, making the split of the 
total energy transfer term into its individual components redundant in this wavenumber 
range. In sec. 7.3.1 we concentrate on the contributions of the different interactions to 
transfers in the inertial range of total energy and set /? = 3, thus taking into account 
only the region in wavenumber space where this scaling is established. Since the values 
of P and P' may influence the direction of energy transfer, a similar approach may be 
useful to calculate energy and helicity transfer at the very low wavenumbers. However, 
this awaits consensus on the low-wavenumber scaling of the magnetic and kinetic energy 
spectra. Furthermore, the integrals must be cut off at some wavenumber such that a 
single scaling exponent for the wavenumber range of interest can be studied. As the 
extent of the inertial range will grow with increasing Reynolds number, contributions 
from the production and dissipation ranges can safely be neglected, as they will become 
very small compared to the extent of the inertial range. However, in the low wavenumber 
region, this argument is not applicable and further work is necessary in order to establish 
if very nonlocal interactions contribute significantly to the transfers of magnetic energy 
and helicity in the low wavenumber range or not. 

7.2. Magnetic helicity transfer 

Using the decomposition into helical modes, the transfer term in the evolution equation 
of the magnetic helicity can be expressed through the transfer term in the evolution 
equation of the magnetic energy, that is 

T«(A:,p,g) = |T«^(fc,p,g) , (7.18) 

and only the transfer term which originates from the induction equation is present, since 
Hmag Is a purely magnetic quantity and as such only implicitly depends on the evolution 
of the velocity field. 

Since the magnetic helicity is an ideal invariant, the transfer term in the spectral 
evolution equation of the magnetic helicity vanishes if integrated over all space, therefore 
similar to the flux of total energy, the flux of magnetic helicity through wavenumber k 
due to a given interaction (f), 

pk poo poo 

= -fdk'J^ T^l{k',p,q)dpdq, 


(7.19) 
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can be written as the sum of two contributions 

1 Q fk 

nH\k)=- T^igik',P,q)dp dq 

1 pk O P'^ /»oo 

~ 2 Jo Jk 

Following the approach explained in sec. 7.1 the integral becomes independent of k and 
one obtains the following expression for the flux of magnetic helicity through k 


nnik) =k‘^-d' dv 



SvTjilgivA,w) 


vd'-^ - 1 
j3' -2 


■ (7.21) 


7.3. Cascades and wavenumber-dependent transfers of total energy and magnetic helicity 

From the expressions (7.17) and (7.21) for the fluxes of total energy and magnetic 
helicity, respectively, it is now possible determine the sign of the fluxes and hence the 
direction of energy and magnetic helicity transfers using the results from the stability 
analysis. If the total energy flux is positive, energy is transferred from smaller to 
larger wavenumbers and if it is negative, energy is transferred from larger to smaller 
wavenumbers. As the magnetic helicity is not positive definite, the situation is slightly 
different. For positive magnetic helicity a positive flux indicates forward transfer just 
as for the total energy. For negative magnetic helicity a negative flux indicates forward 
transfer while a positive flux indicates inverse transfer. However, as this situation is 
symmetric we assume positive helicity throughout the analysis. 

In sec. 6 unstable solutions of (3.3) and (3.4) were interpreted as leading to energy 
transfer out of the unstable mode into the two modes it interacts with for a given helical 
mode interaction (i). If and Bg^ are the unstable modes, this interpretation leads to 

dt\Bg^\^=Tgi^{p,k,q)<Q, (7.22) 

and 

dtWsf = Ti;^(p,fc,g) +T«(p,fc,g) < 0 . (7.23) 

The instability assumption therefore attributes signs to the transfer terms, which will 
determine their respective contributions to the overall energy (and magnetic helicity) 
transfer. Note that dt\UsPf and dt\BsPf cannot have different signs, as both signs are 
determined from the existence of exponentially growing solutions of the system (4.4). 

We now treat the three helicity combinations which lead to unstable solutions sep¬ 
arately assuming without loss of generality that Sp = 1. Having determined the signs 
of the transfer terms within our framework, we now use these results to calculate the 
contributions of the individual transfer terms to the fluxes of total energy and magnetic 
helicity though a given wavenumber. 

7.3.1. Total energy cascades 

For the (inertial range) energy cascade the flux is wavenumber-independent leading to 
/3 = 3 in (7.17). Hence the integrand in (7.17), which determines the sign of the total 
energy flux, becomes 


Ie = (w, 1, u) In w -I- (v, 1, w) In r; , 


(7.24) 
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where we remind the reader of the wavenumber ordering v ^ 1 ^ w. That is, the term 
T(d 

{w, 1, v) describes energy transfer in and out of the largest wavenumber modes while 

T(d 

{v, l,w) describes energy transfer in and out of the smallest wavenumber modes. 

Using the signs of the transfer terms determined for the three helicity combinations 
depending on wavenumber ordering, we can now deduce which helicity combinations 
contribute to forward or inverse cascades of total energy. 

• Sy - - Syy 

For this case we can deduce from the results of the stability analysis summarised in 
tables 1-3 that T(*)(1 ,u, w) < 0, as modes corresponding to the middle wavenumber are 
unstable, while (u, 1, rc) > 0, as modes corresponding to the smallest wavenumber 
are stable and hence these modes can only receive energy from the modes at higher 
wavenumbers. The sign of depends on the values of cross-helicity and the 

ratio of magnetic to kinetic energy. For a magnetically dominated system 1, v) < 0 

and we obtain in this case an inverse cascade of total energy, as 

/£. = (w, l,v)lnw + (u, 1, w) In u < 0 . (7.25) 

If the kinetic energy is much larger than the magnetic energy, cancellations between the 
two terms in Ie occur. The term (w, 1, u) In w is now positive, since the modes at 
the largest wavenumber can only receive energy, thus contributing to a forward cascade. 
For intermediate cases the value of the cross-helicity becomes decisive as high cross- 
helicity quenches the inverse transfer in this case. In summary, we expect inverse cascade 
contributions from this combination of helicities if the magnetic energy dominates, 
while for larger kinetic energy high values of cross-helicity quench the inverse transfer 
contribution to some extent. 

• Sy = Sl ^ Sy] 

From tables 1-3 the instability assumption imposes < 0 and < 

0 as modes corresponding to the middle and largest wavenumbers are unstable. The 
stability of the lowest wavenumber modes now depends on several parameters. However, 
the detailed conservation property 

T^^\w,l,v) + T^^\l,v,w) + T^^'>{v,l,w) = 0 (7.26) 

asserts that if two of the transfer terms are of like sign, then the third one must be of 
opposite sign. Since there are no constraints on the instability of the modes corresponding 
to the two largest wavenumbers, we take 1, w) to be positive and conclude that 

this combination of helicities leads to an inverse energy cascade as Ie <0, and we note 
that this case behaves very differently from its hydrodynamic analogue, where it led to a 
forward cascade of kinetic energy (Waleffe 1992). We also note that this inverse cascade 
should always be present, as it is not subject to constraints from Hc{p) and \UsJ/\Bsj,\. 

• 

Analogously, we obtain < 0, since the modes corresponding to the middle 

wavenumber are unstable. As the stability of the two other transfer terms depends on 
several constraints, no clear assessment is possible. If we assume them to be of like sign 
and thus positive, as not all transfer term can have the same sign due to the detailed 
conservation property, we obtain a contribution towards an inverse cascade. However, if 
we assume them being of opposite sign, contributions to inverse and direct cascades are 
possible. We note that the instability leading to forward transfer in this case is damped 
by high values of Hc{p). 
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7.3.2. Magnetic helicity transfer in the inertial range of total energy 
For /3' = 3, the integrand In in (7.21) becomes 

Ih = T^lg{wA.v)s^(w-1)+ T^lg{v,l.w)sy{v-1) . (7.27) 

Using the signs of the transfer terms determined for the three helicity combinations, we 
can now deduce which helicity combinations contribute to a forward or inverse cascade of 
magnetic helicity. As can be seen in (7.27), there is an explicit dependence of the magnetic 
helicity flux on the helicities of the interacting modes. In the following we assume si = 1. 

• S'D = Si = Syj 

The integrand Ih becomes 

lH = T^ig{w,l,v){w-l)+Tjfl{v,l,w){v-l) . (7.28) 

As the signs of the magnetic energy transfer term deduced from the stability analysis 
are the same as for the total energy and In w and w — 1 are both positive while In v and 
u — 1 are both negative, the result for the helicity transfer reflects the results for the total 
energy cascade, thus for this helicity combination total energy and magnetic helicity will 
be transferred in the same direction, which can be both forward and inverse in this case. 

• Sy = Sl ^ Syy 

The integrand Ih becomes 

lH = -T^Uw,l,v){w-l)+T(fl{v,l,w){v-l) , (7.29) 

where the contributions from the largest wavenumber modes now enter with the opposite 
sign. Compared to the total energy flux, which was purely inverse in this case, we obtain 
the possibility of simultaneously a forward helicity flux and an inverse energy flux. 

• ^ ^1 - 

The integrand Ih becomes 

lH = T^ig{w,l,v){w-l)-T^lg{v,l,w){v-l) , (7.30) 

where the contributions from the smallest wavenumber modes now enter with the 
opposite sign. Compared to the total energy cascade, again we find that it is possible 
to have a transfer of magnetic helicity in the opposite direction to the transfer of total 
energy. 

In this subsection we determined the direction of the magnetic helicity transfer in the 
inertial range of total energy for different combinations of helicities and compared the 
results to those for the total energy cascade. We found that a cascade of total energy 
is possible in one direction while the transfer of magnetic helicity may proceed in the 
opposite direction. 

7.3.3. Magnetic helicity cascades 

In the inertial range of magnetic helicity the flux of magnetic helicity is wavenumber- 
independent resulting in /?' = 2 in (7.21). Therefore the integrand Ih in (7.21) becomes 

Ih = r^*lg(w, I,u)s™lnu;-hr^*ig(u,I,w)s„lnu . (7.31) 

For the three different helicity combinations this leads to 

• Sy Syy 

The integrand in this case is of the same form as the integrand Ie for the total energy 
cascade (that is, if /? = 3 in I^;) 

Ih = T^i„{w,l,v)\nw+ {v,1,w)\tlv , 


(7.32) 
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hence the results for the cascades of magnetic helicity are the same as for the cascasdes 
of total energy. 

• Sy = Sl ^ Syy 

The integrand in this case has a different form compared to the integrand Ie for the 
total energy 

Ih = -Tjillg{w,l,v)lnw+ T^lg{v,l,w)\nv , (7.33) 

hence the results for the cascades of magnetic helicity differ from the total energy 
cascades. In particular, this case may lead to a nonhelical reverse energy transfer while 
the helicity cascade may be forwards, due to the contribution from Tmag(w,l,v) now 
having the opposite sign in Ih compared to Ie- 

* ^ ^1 - 

Again, the integrand in this case has a different form compared to the integrand Ie for 
the total energy 

Ih = Tj^ig{w,l,v)\nw -T^lg{v,l,w)\nv , (7.34) 

hence the results for the cascades of magnetic helicity are different from the total energy 
cascades. In particular, this case may lead to a nonhelical reverse energy transfer while 
the helicity cascade may be forwards, due to the contribution from Tmag(v,l,w) now 
having the opposite sign in Ih compared to Ie- 

7.3.4. Magnetic energy transfer in the inertial range of magnetic helicity 

For j3' = 2, the contributions to the integrand Ie due to magnetic energy transfer are 



The signs of T^ag and T are the same by eqs. (7.22)-(7.23), and In re and w — 1 are both 
positive while Inu and v — 1 are both negative. Hence, the result for the contributions 
of these terms to the total energy transfer in the inertial range of magnetic helicity 
is the same as in the inertial range of total energy for all helicity combinations. That 
is, magnetic energy transfer and conversion in the inertial ranges of total energy and 
magnetic helicity proceed in the same direction. 

This assessment of contributions to forward and inverse transfers and cascades is 
based on an analysis of the nonlinear terms in the MHD equations only, thus neglecting 
the symmetry-breaking effect of dissipation creating an energy sink at the small scales. 
Accounting for this effect, it is plausible that the contributions from transfer terms leading 
to forward transfer are higher weighted than contributions leading to inverse transfer. 
This is particularly relevant in interactions where forwards and reverse contributions 
are present and the overall transfer depends on cancellations between the two terms. It 
would perhaps be safest to attribute these cases to forwards rather than inverse energy 
cascades. 

Although it is not possible to exactly determine which helical interactions are higher 
weighted than others, some information can be obtained from the magnitude of the 
geometric factor gkpq defined in eq. (3.5). The magnitude of gkpq depends on the helicity 
combinations since it involves the helicity-dependent factor I = Skk + SpP+Sqq. Therefore 
it parametrises the strength of a given helical interaction, and the case of all helicities 
being of the same sign gives the largest value of |/|, since in this case |/| = \k + p + q\- 

For the reverse transfers, that is, for k < p,q, the factor |/| takes the smallest value 
for the case Sk = Sp Sq, since |/| = \k + {p — g)]. Note that in this case I becomes 
small for small k even in the nonlocal limit k « p q, suggesting that the nonhelical 
reverse transfer found in this case is less efficient in increasing spectral power at the 
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very low wavenumbers. The remaining class of helical interactions Sk Sp = Sq leads to 
|/| = \k — {p + q)\. In this case |/| does not necessarily become small for small k which 
is due to the contribution of nonlocal interactions, where p and q are large compared 
to k. According to the results from the stability analysis, in the nonlocal limit unstable 
solutions occur for the case Sk = Sp = Sq only if \Bs^ \ > |C4p| and for the case Sk ^ Sp = 
Sqii\UsJ>\BsJ. 

It is therefore possible to deduce within the framework of the instability assumption 
that most of the increase in energy at the very largest scales (in a magnetically dominated 
system) is mainly due to a breaking of mirror-symmetry, which had been established 
before by Frisch et al. (1975) using a different approach. That is, it is due to the presence 
of kinetic and magnetic helicity, since interactions of the type Sk = Sp = Sq, which account 
for most of the inverse transfer, can only occur in significant numbers for fields consisting 
of many modes with the same helicity. Recent numerical results in hydrodynamics showed 
that there is an overall reverse flux of energy only when the system mainly contains helical 
modes of the same sign. As soon as a small amount of oppositely-polarised modes are 
introduced, the usual direct cascade is recovered (Sahoo et al. 2015). 

In summary, in this section we determined the direction of total energy and magnetic 
helicity transfers in their respective inertial ranges. Not surprisingly, we found that fully 
helical magnetic fields lead to inverse cascades of magnetic helicity and magnetic energy, 
but the analysis also showed that an inverse energy cascade is possible for nonhelical 
magnetic fields, which is a new theoretical result. However, due to the coupling of the 
momentum and induction equations, within this framework it is not possible to determine 
the nature of the energy transfers resulting from an instability of a given steady solution, 
since the same eigenvalue controls the growth of the exponential solution of (4.4) for both 
the magnetic and the velocity held. Nevertheless, for some special cases the evolution 
equations (3.3)-(3.4) decouple and more detailed information becomes available. These 
cases are treated in the following section. 

8. Special solutions and the (kinematic) dynamo 

Having established the general case, we now draw attention to special cases where 
the analysis becomes much simpler and which are relevant to specihc problems in MHD 
such as the kinematic dynamo. In sec. 4.2 we analysed the stability of general steady 
solutions of the dynamical system (3.3)-(3.4), which describes the evolution of a triad of 
interacting helical modes. Using the notation (4.1), the general steady solutions were of 
the form (0, Usj,,0] 0, 0). In this section we now study the cases where either Ug^, = 0 

or Hsp = 0, that is we analyse the stability of steady solutions of (3.3)-(3.4) of the form 
(0, C/gp, 0; 0,0, 0) and (0,0,0; 0, Bsp,0). The former case may be of particular interest due 
to its relation to dynamo action. 

8.1. The kinematic dynamo 

For small magnetic fields the Lorentz force is small compared to inertial forces, and 
can be neglected in the momentum equation. This decouples the momentum equation 
from the induction equation and defines the kinematic dynamo problem. In our setting, 
it corresponds to \Us^\/\Bg^ \ >> 1, and terms proportional to |Rgp| can be neglected as 
they are very small compared to terms proportional to | Ug^ \. 

Alternatively, one could also consider the steady solution Bg^ = 0 while Ug^ ^ 0, 
which would correspond to a stability analysis of a flow held at a particular length scale 
subject to small perturbations of the magnetic and velocity helds, where the magnetic 
Held perturbation may be viewed as the magnetic seed Held to be ampHHed by dynamo 
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action. In this setting we observe from (3.3) that the term corresponding to the Lorentz 
force disappears while in (4.3) terms involving Bg disappear, thus the system simplifies 
to 

9tUs^ = \gkpq\'^{spP - Sgq){skk - Spp) {UgJ'^Us^ , ( 8 . 1 ) 

dtbs^ = -\gkpq\^Skksqq\UgJ'^bg^ . ( 8 . 2 ) 

As the only contribution to the evolution of the magnetic field now comes from the 
velocity field, we associate the remaining terms in (3.4) with dynamo action. From (8.2) 
we observe that this system has exponential solutions leading to magnetic field growth 
if Sk ^ Sq, regardless of wavenumber ordering. So for energy transfer from Ug into 5^^ 
(and bg^) to become possible, the magnetic modes at wavenumbers k and q should be of 
opposite helicity. 

For small k, nonlocal interactions with k « p q provide most transfer into bg^. This 
is because the eigenvalue determining the growth of the exponential solution of (8.2) is 
larger for q » k than for q k, thus the perturbations should grow faster in the former 
than in the latter case. Hence, according to the instability assumption, Ug^ loses energy 
in favour of bg^, mainly due to nonlocal interactions if 5^^, describes the largest scales of 
the system. 

8.1.1. The a-effect 

One well-known example of a large-scale dynamo is the a-effect of mean-field electro¬ 
dynamics (see e. g. Moffatt (1978)), where a is a coefficent in the mean-field induction 
equation related to the kinetic helicity of the flow. The a-effect leads to a generation 
of large and small-scale magnetic helicities of opposite sign (Brandenburg 2001, 2003). 
A positively helical velocity field generates magnetic field perturbations leading to the 
large-scale component of the magnetic field becoming negatively helical. By conservation 
of magnetic helicity, the small-scale component of the magnetic field then has to become 
positively helical (and more so if the initial magnetic field was positively helical). That 
is, the small-scale magnetic and kinetic helicities are of the same sign. 

It is plausible that the type of interaction {0,Ug^,0;0, Bg^ = 0,0) for k < p,q with 
Sk 7 ^ Sp = Sq can be associated with the a-effect. First, nonzero small-scale kinetic 
helicity (we have Sp = Sq) is present. Second, the magnetic field growth at the large 
scales is described by (8.2), where magnetic fluctuations at k and q of opposite helicities 
are necessary to obtain an unstable solution. That is, the large-scale magnetic field has 
opposite helicity to the small scale one, reminiscent of the a-dynamo. We also note 
that this combination of helicities produces a transfer of kinetic energy from small to 
large scales (Waleffe 1992). Thus this type of interaction feeds into the magnetic and 
velocity fields on scales larger than the characteristic scale L = 1/p of the velocity field. 
The magnetic field mode which is amplified by this process has helicity opposite to the 
velocity field at p, which conforms to expectations in terms of the a-effect. 

8.2. Excitation of a flow by the Lorentz force 

For the other special solution (0, 0,0; 0, Bg^ 0,0) the system (4.2) - (4.3) simplifies 
to 

= \gkpq\ {SpP — Sqq)Sqq\Bgfl Ug^ , (8-3) 

~ Skki^Skk Spp) bg^ , (8-4) 

and we note that the inertial term in (3.3) and the ‘dynamo’ term in (3.4) are now 
absent and the system of coupled ODEs has split into two decoupled ODEs. This case 
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may perhaps be associated with the generation of turbulence caused by the action of the 
Lorentz force on the fluid (i.e. energy conversion from to Ug,. or Ug^) and interscale 
transfer of magnetic energy from Bg into or bg . Exponentially growing solutions of 
(8.4) only occur if Sp = Sk and k < p, leading to a reverse transfer of magnetic energy. 
Exponentially growing solutions of (8.3) occur ior p > q and Sp = Sq leading to forward 
and reverse transfers corresponding to k > p and k < p, respectively. Interestingly, 
energy transfer only becomes possible if the magnetic field is helical and the helicity of 
the velocity field mode does not affect the analysis. 


9. Conclusions 

The four main results of the present work are: First, unlike in non-conducting fluids 
(Waleffe 1992), the stability analysis shows that in MHD turbulence energy can be 
transferred away from the smallest scales in a triad interaction. Second, the stability 
analysis reveals mechanisms for reverse energy transfer for nonhelical magnetic fields, 
in which case it does not need to be driven by the inverse transfer of magnetic helicity. 
Third, forward energy transfers are more quenched in regions of high cross-helicity than 
reverse energy transfers. Fourth, we expect significant cancellations to occur between the 
contributions to forward and reverse transfers, as on several occasions they occur with 
opposite signs in the same equation. Our theoretical analysis was conducted within the 
framework of the instability assumption, and it is crucial to discuss the results within 
the wider context of MHD turbulence research. 

Interscale energy transfers between the two different vector fields as well as within the 
same fields have been studied by several groups for freely decaying (Debliquy et al. 2005; 
Brandenburg et al. 2015) and stationary (Brandenburg 2001; Alexakis et al. 2005; Carati 
et al. 2006; Cho 2010) MHD turbulence as well as for the kinematic dynamo regime 
(Mininni et al. 2005) and for magnetic helicity transfer (Alexakis et al. 2006), using 
shell-filtered transfer terms calculated from DNSs or from a helical shell model (Stepanov 
et al. 2015). In the stationary case, it was found that transfers between the same fields 
are mainly local while transfers between different fields were nonlocal, and transfers from 
the injection scale to the largest scales in the system were observed. In the decaying case, 
energy transfers were generally found to be mainly local. However, transfers between 
different fields were more nonlocal than transfers between the same fields. Furthermore, 
large cancellations occurred between the contributions to forward and reverse transfers 
(Debliquy et al. 2005). The analysis presented here also predicts cancellations between 
these contributions to occur, thus being consistent with the aforementioned numerical 
results. 

In terms of locality and nonlocality of energy (and helicity) transfer, we found that 
nonlocal interactions contribute to forward transfer only for interactions of helical modes 
with unlike helicity and mainly if the kinetic energy exceeds the magnetic energy. 
Interestingly, for inverse transfers we find less constraints on nonlocal interactions. In 
particular for magnetically dominated systems nonlocal interaction between modes of 
like helicity contribute to reverse energy transfer. In view of the cancellations that occur 
between forward and reverse transfers, the inverse cascade may thus have a significant 
nonlocal component which is not cancelled by forward transfers within the same triad 
interaction. 

A numerical study of large-scale magnetic field generation in helically forced isotropic 
MHD turbulence was carried out by Brandenburg (2001). It was found that the injection 
of energy from the velocity field into the magnetic field occurs directly from the forcing 
scale into the largest resolved scale, implying that this is a nonlocal process. Due to the 
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non-locality of the observed increase in spectral power of the magnetic field at the lowest 
resolved wavenumber k = 1 and the excellent agreement of numerical results with an 
a-dynamo model, the transfer of energy into the k = 1 mode is explained by the a-effect 
rather than an inverse cascade, and it is shown to occur after saturation of the small-scale 
dynamo. Our results in sec. 8 suggest that one type of helical mode interaction may be 
mapped to the a effect, and we established that large-scale dynamo action is more active 
in the nonlocal limit than certain other types of interactions. 

One of the main results of the present work is the possibility of inverse energy transfer 
for nonhelical magnetic fields. Such inverse transfer has recently been found in high 
resolution DNSs of slightly compressible (Brandenburg et al. 2015) and relativistic (Zrake 
2014) MHD turbulence. An analysis of interscale transfers showed that this inverse 
transfer was mainly due to energy transfer away from the medium scale (see Supplemental 
Material of Brandenburg et al. (2015), last figure), while energy transfer away from the 
smallest scales also occurred. The analytic approach put forward in the present paper also 
shows that energy is transferred away from the medium and small scales for interactions 
of modes with unlike helicities, thus being qualitatively consistent with these numerical 
results. However, since no numerical work decomposing the MHD equations into helical 
contributions as suggested by Biferale et al. (2012) and Biferale & Titi (2013) has been 
carried out so far, we are not in a position to claim numerical confirmation of our results. 

Having discussed numerical results, we now turn to measurements of energy transfer 
in the solar wind. Unlike in our own analysis and in the numerical results discussed 
so far except for (Cho 2010), a background magnetic field is present in the solar wind. 
Recent measurements at 1 AU (Stawarz et al. 2010) showed negative Elsasser fluxes in 
regions of high cross-helicity, giving possible evidence of inverse energy transfer in these 
regions, which cannot be explained by selective decay as cross-helicity cascades forwards 
(Frisch et al. 1975). Our results may be helpful in explaining this phenomenon as one 
of the results we obtained was a quenching of forward energy transfer in regions of high 
cross-helicity, leaving more inverse transfer to perhaps dominate the dynamics in these 
regions. 

In subsequent work (Coburn et al. 2014) concerns were raised on the implications of the 
effect of expansion in the solar wind especially in regions of high He. Expansion effects 
had been neglected in the previous analysis. The authors restrict their measurements 
to regions where the relative cross helicity CTc is not too large, that is 0 < |(Jc| < 0.5 
and measure positive energy fluxes on average, while the instantaneous flux shows large 
variations including negative values. It is shown that the broad distribution of the 
measured instantaneous fluxes are related to intermittency of the energy cascade in terms 
of the variability of the energy flux (Politano & Pouquet 1995; Karimabadi et al. 2013) 
and not caused by experimental uncertainty. The various possibilities of energy transfer 
in forward and reverse directions determined in the present work are consistent with these 
measurements, as they also would result in broader tails of the probability distribution 
of the energy flux, even if on average energy transfer proceeds in the forward direction. 
As for the concerns about the validity of the negative fluxes measured by Stawarz et al. 
(2010), our results do suggest that the measured inverse fluxes may be a genuine effect 
due to quenching of forward energy transfers if He is large. 

Since most of this discussion is based on statements of plausibility rather than certainty, 
more work clearly has to be carried out before a decisive result on energy transfer in MHD 
turbulence can be achieved, and we hope that our analysis constitues a step forwards in 
this direction. As suggested by Biferale et al. (2012); Biferale & Titi (2013), one could 
study energy and helicity transfers numerically by projecting out helical modes of a 
particular sign, similar to work done by Biferale et al. (2012); Biferale & Titi (2013) 
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and Sahoo et al. (2015) in hydrodynamic turbulence. However, numerical verification 
of reverse spectral transfer due to the particular nonhelical interactions found in the 
present work may be difficult to obtain in that framework, and a particular DNS study 
concentrating on inverse transfer for nonhelical magnetic fields using the full MHD 
equations subject to small-scale forcing may be needed in order to provide further insight. 
An analysis of Fourier-filtered transfer terms from DNSs of highly unbalanced MHD 
turbulence compared to balanced MHD turbulence could be carried out in order to verify 
(or otherwise) the proposed quenching of forward transfers by high values of He, especially 
as it is not possible to quantify this effect from theoretical analysis only. On the analytical 
side, the present work may be extended to include the effects of a background magnetic 
field and of compressive fluctuations, which would be included in the decomposition 
of the velocity field as modes parallel to the wavevector k. Asides from providing an 
advance in MHD turbulence research on a fundamental level, this work may contribute 
to the further theoretical understanding of various physical processes involving inverse 
cascade and dynamo action. This includes the evolution and origin of cosmological and 
galactic magnetic fields as well as solar physics and the dynamics of laboratory plasmas 
and turbulence in liquid metal flows. 
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Appendix A. The dependence of Q on the cross-helicity 

In section 4, the parameter Q was defined as Q = a6 — (dj, where a, j5, 7 and 6 were 
the entries of the matrix in (4.4). Using the expressions for these terms given in (4.8), 
we obtain 

Q= \9kpq\'^SkkSqq{Skk - Spp){Sqq- Spp) [\UsX + \BsX ■ (^1) 

However, in (4.12), instead of the term Re([[/*pi?sp]^) a term involving Hc{p) appeared. 

In general, the helical coefficients and are related by a complex number M = 
m+in such that Bg^ = MUg^. Expressions for m and n can be found by decomposing the 
two fields into their real and imaginary parts. Let Ug = Ui + iU 2 and Bg = Bi + ii? 2 - 
Then 

m =—^(UiHi-b [/ 2 R 2 ) and (A2) 

n=Tf^iUiB2-U2Bi), (A3) 

and we note the constraint = \Bg P/|Us p — which follows from the definition of 
M. Decomposing the cross-helicity in the same way results in Hc{p) = \Ug pm. Now we 
can relate Re([[/*^i3sp]^) to the cross-helicity by rewriting it in terms of the components 



of Usj, and Bs/. 
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=(C/iRi + U^B^f - {U1B2 - U^B^f (A 4 ) 

= \UsXm^-\UsX (A5) 

=2HM^-\U,f\B,f , (A6) 

and we obtain (4.12) by substitution of this expression for Re([C/*^i3sp]^) into (Al) 
Since the maximum and minimum values of \H(.{p)\ are |C7sp||5sp| and 0 respectively, 
it is useful to define the relative cross-helicity p = Hdp)/{\Us^\\Bs^\), which takes values 
between -1 and 1. We obtain 

Re(KpS«.]") = \UsX\B,^\\2p^ - 1) , (A7) 

which is bounded by — | b"sp P 1 73sp P and | t^sp P | Bsp P, where the first value is the case of 
vanishing cross-helicity and the latter occurs when there is maximal cross-helicity. This 
implies that the term + \BsJ^ + 2\UsJ'^\BsJ‘^ — AHdp)^) in (4.12) cannot be 

negative. 


Appendix B. — Q > 0 for sp = Sq ^ Sp 

In section 5.1, the result was dependent on whether — Q is positive or negative. 
Recall that the helicity combination in question was Sk = Sq ^ Sp. Therefore — Q 
becomes 

x^ -Q =^[\UsX{k<l +{k+p){q+p)f + \B{k(k + p) + q{q + p)f 
+ 2\UsX\BsX\kq+{k+p){q+p)\[k{k + p)+q{q+p)\] 

- [kq[k+p){q+p)]{\UX^ + \BsX - 2 Re([C/;i 3 .p] 2 )) 

= \[\UsXikq-{k+p){q+p)f -f \BsX{k{k + p) - q{q + p)f 
+ 2\UsX\BsX[kq+{k+p){q+p)\[k{k + p)+q{q+p)\] 

+ 2Re{[U:^BsX)[kq{k+p){q+p)] ■ (B 1) 

If Re([[/*pi?sp]^) > 0, then this implies — Q > 0. If Re([I7*^Rsp]^) < 0, then some 
more steps are required to show that — Q > 0. In general, we have |Re([I7*^Rsp]^)| ^ 
\UsX\BsX appendix A), hence assume Re([I7*^Rsp]^) = —\UsX\BsX this would 
be the most negative value this term can take. It corresponds to zero cross-helicity at p. 
Equation (B 1) is now an inequality and reads 

x'^ [\UsXikq- {k + p){q + pX + \BsXik{k + p) - q{q + pX] 

+ ^\UsX\BsX[kq + (k + p){q + p)][kik + p) + q{q + p)] 

-2\UsX\BsX{kqik+p)iq+p)] ■ (B2) 

The result a;^ — Q > 0 follows immediately if we can show 

[kg +{k+p){q+ p)] [kjk + p) + gjq + p)] 
kq{k + p){q + p) 


(B3) 
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Defining v = p/k and w = q/k^ this inequality can also be written as 


1 

V + w 


w 1 + u + w(v + w) 

^^^ > 4 
1 + ri w 


(B4) 


where v and w must satisfy l^u^w<l+v due to the triad geometry. This implies 
that the first and second terms on the left-hand side (LHS) of the inequality are smaller 
than unity, and the third term is the largest of the three. In order to prove the inequality, 
we therefore concentrate on the third term and derive a lower bound for it. We proceed 
by defining a family of curves Fy(w) by 


Fy{w) 


1 -b V -b w{v + w) 

W 


(B5) 


and aim to find their minima depending on the parameter v. The triad geometry will 
then give a smallest allowed value of v and thus provide the smallest possible value of 
the third term on the LHS of (B4). Differentiating Fy{w) with respect to w yields 


dw 


Fy{w) 


w{2w -b u) — (1 -b U -b w[v + w]) 


— (1 -b v) 




(B6) 


setting this expression equal to zero results in w = ±-\/l -b v, and it is easily checked that 
{d^/dw'^)Fy{w) > 0. By definition w > 0, hence w = \/l + v, and we obtain 


7^ / n -X 1 -b r’ -b -\/l -b v{v + \/l + v) r. - 

Fy{\/l+v) = - , ^ ^-- = 2VI+V + V . (B7) 

Vl + v 

This expression is minimal for ri = 1 and in this case equal to Fi(-\/2) = 2^/2 + ! ~ 3.828, 
which is the lower bound for the largest term on the LHS of (B 4). Substituting the 
corresponding values v = 1 and w = \/2 into the remaining terms on the LHS of (B 4) 
we obtain 

—-by2-b2v^-b 1 ~ 1.828-b 3.828 > 4 , (B8) 

1 + ^2 

thus the inequality (B 3) is satisfied and — Q > 0. 


Appendix C. Graphical determination of contraints on stability 

As explained in the main body of the text, the term Q given in (4.12) determines 
the stability of the system (4.4) if a; < 0. As such, a solution is unstable if Q < 0 or 
a — Q < 0, where the latter case is the more difficult to determine, as the sign of 
x^ — Q depends on the shape of the wavenumber triad, the cross-helicity and the ratio 
\Usp\/\Bsj,\. Given the multitude of possibilities that can emerge for this, the simplest 
way of determining the constraints on the stability of a solution of (4.4) is using a 
graphical method. For each combination of helicities we plot x^ — Q for several set values 
of \Usp\/\Bsp I and F[c(p) in order to show in which parameter range instabilities are more 
likely to occur. 

The dependence oi x^ — Q on the triad k,p, q can be reduced to a dependence on the 
triad’s shape by rescaling each wavenumber similar to the procedure in appendix B, which 
enables us to use two-dimensional plots and the triad geometry to obtain the necessary 
information. Figures 1-4 show the function x^(v,w) — Q(v,w) for the different cases 
shown in tables 1-3, where v and w correspond to the smallest and largest wavenumber 
in a given triad, rescaled by the middle one such that the triad geometry enforces the 
contraint 0<ri^l^r(;<l-bw, hence each wavenumber pair {v, w) describes a shape of 
triad. Each subfigure corresponds to set values of F[y{p) and \Us^\/\Bs^\, while each point 
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Figure 1: Plots of f{v, w) = — Q for various values of \UsJ/\Bsp \ and cross-helicity for 

case 1 in appendix C (s^ ^ Sp = Sq,k < p < q). The upper grey triangle is ruled out by 
the condition w < 1 + n and unstable values are shown in white. The ratio |{7^^ |/1| 
increases from left to right, with each column of subfigures taking the values 0.01, 0.1 
and 1 respectively, while each row takes the following values of relativ cross-helicity: 
^c(p)/(|C/.J|i3.J) = 0,0.5, 0.9 and 1. 


(v, w) in a particular graph corresponds to a class of triad interactions characterised by 
their shape. Regions in wavenumber space excluded from the analysis by the constraints 
of the triad geometry are shaded in grey, positive values of — Q leading to stability are 
indicated in black and negative values oi — Q leading to unstable solutions are marked 
white. Across the main four figures, Hc{p) increases towards the bottom of the figure 
while |C 7 sp |/|.Bsp I increases from left to right, leading to the constraints summarised in 
tables 1-3. Depending on the wavenumber ordering, the definitions of v and w are slightly 
different, and we explain the procedures for each case individually. 

• Sk ^ Sp = Sq and k < p < q 

In this case we rescale all wavenumbers by p, such that v = k/p and w = q/p. As can 
be seen in fig. 1 , for decreasing \Usp\/\Bs^\ and increasing Hc{p) less and less unstable 
solutions occur and we obtain the constraints on split transfer shown in table 1 . 

• Sk ^ Sp = Sq and p < k < q 
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Figure 2: Plots of f{v, w) = — Q for various values of \UsJ/\Bsp \ and cross-helicity for 

case 2 in appendix C (s^ ^ Sp = Sq,p < k < q). The upper grey triangle is ruled out by 
the condition w < 1 + n and unstable values are shown in white. The ratio |{7^^ |/1| 
increases from left to right, with each column of subfigures taking the values 0.01, 0.1 
and 1 respectively, while each row takes the following values of relative cross-helicity: 

= 0,0.5, 0.9 and 1. 


In this case we rescale all wavenumbers by fc, such that v = p/k and w = q/k. As can 
be seen in fig. 2, for decreasing \Usp\/\Bs^\ and increasing Hc{p) less and less unstable 
solutions occur and we obtain the constraints on forward transfer shown in table 3. 

• Sk = Sp ^ Sq and p < k < q 

In this case we rescale all wavenumbers by k, such that v = p/k and w = q/k. As can 
be seen in fig. 3, for decreasing \Us^\/\Bg^\ and increasing Hc{p) less and less unstable 
solutions occur and we obtain the constraints on forward transfer shown in table 3. 

• Sk = Sp = Sq and k < q < p 

In this case we rescale all wavenumbers by g, such that v = k/q and w = p/q. As can be 
seen in fig. 4, now for increasing \Usj,\/\Bsp \ and increasing Hdp) less and less triads lead 
to unstable solutions and we obtain the constraints shown in table 2 on reverse transfer. 
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Figure 3: Plots of f{v, w) = — Q for various values of \UsJ/\Bsp \ and cross-helicity for 

case 3 in appendix C (s^ = Sp ^ Sq,p < k < q). The upper grey triangle is ruled out by 
the condition w < 1 + n and unstable values are shown in white. The ratio |{7^^ |/1| 
increases from left to right, with each column of subfigures taking the values 0.01, 0.1 
and 1 respectively, while each row takes the following values of relative cross-helicity: 

= 0,0.5, 0.9 and 1. 


Appendix D. Similarity scaling 

In the respective inertial ranges 

Ekiniak) / Ekin{k) = a“”, E^ag{ak) / E^ag{k) = a”™ (D 1) 

where a is a real number and n > 0 and m > 0 are the spectral indices of the kinetic 
and magnetic energy spectra, respectively. From 

EkUk)dk=l [ {\u+{k)\^ + \u_{k)\^)dk , (D2) 

^ J\k\=k 

Emaaik)dk = \ ( (|6+(fc)P + \h^{k)\^)dk , 

^ J\k\=k 


(D3) 
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0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 



0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 



0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 


Figure 4: Plots of f{v,w) = — Q for various values of \Usp\/\Bs^\ and cross-helicity 

for case 4 in appendix C (s^ = Sp = Sq,k < q < p). The upper grey triangle is ruled 
out by the condition n; < 1 + n and unstable values are shown in white. The ratio 
\Us \/\Bs I increases from left to right, with each column of subfigures taking the values 
1, 10 and 100 respectively, while each row takes the following values of relative cross- 
helicity: Hc{p)/{\Usp\\Bs^\) = 0,0.5,0.9 and 1. 

we then find a scaling of the helical coefficients 

Ms(afe) = 6s(afc) = . (D4) 

From (D 1) the scaling of T^jj{k,p, q), T^p{k,p, q) and Tm\g{k,p, q) is then given by 


ap, aq) 

- n--(l+3")/2 - 

(D5) 

T%{k,p,q) 
Tpp(ak, ap, aq) 


(D6) 

Tplik^p, q) 
Tmag{ak, ap, aq) 

Oc \Ac ^ 

(D7) 

T^lg{k,p,q) 

C.X Oc ^ 
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In hydrodynamics n = 5/3, while in MHD there are different predictions for the spectral 
exponent, either m = 3/2 (Iroshnikov-Kraichnan) or m = 5/3 (Kolmogorov). Note that 
n = m = 5/3 implies j5' = j3 = Z while n = 5/3 and m = 3/2 implies jd' = 2 + 5/6. In 
both cases /3 — 2 > 0. 
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